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Abstract. We investigate the space oi N x N dimensional density matrices. We 
show that there exist strata such that boundary states pp with p zero eigenvalues lie 
on or outside the spheres with radii Vp — ^Jp/N{N — p). Moreover, we show that if in 
a certain direction there is a boundary state with q = N — p equal eigenvalues, then 
in the opposite (antipodean) direction exists a boundary state with p = N — q equal 
eigenvalues. 



PACS numbers: 02.10.Yn, 03.65.-w 

Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 

The present renewed interest in the geometry of density matrices describing A^-level 
quantum systems is largely due to the advances in the fundamentals of quantum 
mechanics and in quantum information theory [1] . The recent monograph by Bengtsson 
and Zyczkowski [2] gives an excellent review of the subject and also very extensive 
bibliography. Nevertheless, there are still many important open questions. It is relatively 
easy to account for the required hermiticity and normalization of density matrices. But 
the requirement of positivity (strictly speaking: semi-positivity) causes serious problems. 
One of the methods of investigation of density matrices uses a concept of generalized 
Bloch vectors. This is due to the fact that the space of density matrices is in one- 
to-one correspondence with the space of these vectors (see, for example O HI |5l [6]). 
The recent paper by Kimura and Kossakowski [7j proceeds along these lines. These 
authors construct the spherical-coordinate point of view in the space of Bloch vectors. 
To our minds, perhaps the term directional would be better. Kimura and Kossakowski 
derive inequalities which must be satisfied by the lengths of Bloch vectors in the specific 
directions. Moreover, they show that if there is a pure state in certain direction then in 
the opposite one there occur states with at most one zero eigenvalue. 
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The aim of this work is to simphfy and generahze the results of Kimura and 
Kossakowski. We stress that we deal only with A^-dimensional quantum systems 
described by x density matrices. 

In the next section we briefly present the basic theoretical notions necessary in our 
analysis. We will formulate the expressions giving the metric in the space of density 
operators and hence the tool to measure distances which can also be expressed in terms 
of lengths of Bloch vectors. In the third section we will apply introduced ideas and show 
that the space of density matrices is stratified (that is possesses, colloquially speaking, 
onion- like structure). The fourth section is devoted to a presentation of the directional 
representation of density matrices. This is a brief review of some results of Kimura and 
Kossakowski [7j which seem to be essential for the completeness of this work. In the 
fifth section we discuss and generalize the antipodean properties anticipated in |7]. In 
the last section we give some additional comments and remarks. 

For further convenience, we present here a very simple lemma which is due to 
Harriman [8j and is quite useful in further discussion. 

Lemma. Let Oj E M., {j = 1,2, . . . ,n) and Yl^=i ~ ^ (note that numbers aj need 
not be positive). Denoting J2]=i ^"j = ^) has: 

(i) A>l/n; (la) 

(ii) A = l/n ^ aj = l/n, (j = 1, 2, . . . , n). (lb) 

Proof: Consider the sum S defined as follows 

where we used the condition imposed upon the sum of numbers aj. Then, we have 
A = S + 1/n. Sum S" is a sum of squares so it is nonnegative, then A > 1/n, which 
proves ffToj) . Let us now assume that A = 1/n. Then relation ([2]) entails 5 = 0, which 
means that all aj = 1/n. Conversely, taking oj = 1/n we see that ([2]) gives A = 1/n. 
This completes the proof of the lemma. 

2. Theoretical framework 

Let us denote the maximally mixed state as 

/ 1 1 1 \ 1 . . , 

Pma.=diagi^-,-,...,-j=-l^, (3) 

where Iat is a unit NxN matrix. Hilbert-Schmidt distance between an arbitrary density 
matrix and Pmax is given as 

d'^ {P, Pmax ) = Tr { (p - p^rnax ) ^ } • (4) 

The trace is invariant with respect to unitary transformations and so is pmax- Then 
we choose such a representation in which p = diag(Ai, A2, . . . , Aat), where Xj are 
the eigenvalues of the density matrix and satisfy the requirements: (i) Xj G M 
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(hermiticity), (u) J2j'=i^j ~ ^ (normalization), {Hi) Xj > (semi-positivity). Using 
this representation we have 

j=l \ / j=l 

Since Tr{p^} < 1 we immediately get 

,2/ \ ^ - 1 / ^ 

d [p,pmax) < ^ ■ (6) 

This means that all density matrices lie on or within a sphere (in the space of N x N 
hermitian and normalized matrices) with the center at pmax and radius 



We will call this sphere a large one. On the other hand, it is well known that for > 3, 
not all matrices within this sphere are true density ones. The sphere also contains 
hermitian and normalized but nonpositive matrices. It is also known (see [8]) that there 
exists a small sphere (centered also at Pmax) with a radius 



- Y N{N-iy 

inside which there are only density matrices. The nonpositive ones can be found only 
outside a small sphere. 

Above considerations were made by reference to the eigenvalues of the density 
matrices. Similar ones can be done in terms of the generalized Bloch vectors which are 
introduced as follows. 

The basis in the space of hermitian N x N matrices consists of A^^ elements. We 
take the one consisting of the unit matrix i^v and of A^^ — 1 traceless, hermitian and 
orthonormal matrices {Tj}^^^^. There are various possibilities of the choice of such a 
basis, see the review [6]. We assume that orthogonality holds in the Hilbert-Schmidt 
sense, ie., Ti {TjTk} = 6jk- Then, an arbitrary density operator for a A^-level system 
can be written as 

p = l^^^+ E ^.T, = li^ + v-f, (9) 

i=l 

where (Vi, V2, . . . , Vj\f2_^i) G M^^"^ is called a generalized Bloch vector. Representation 
ensures normalization and hermiticity. Requirement of positivity imposes restrictions 
upon vector V. To our knowledge, their general form is not known although some 
particular cases (SHU [5] were studied in much detail. Discussion of this point, however, 
goes beyond the scope of the present work. 

Bloch vector can be used to express the distance between an arbitrary density 
matrix and the maximally mixed one 

Af2-1 

d\p,Pmax) = Tr {(p - p^axf] = ^ ViVjSij = \V\^ . (10) 

«J=1 
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as it follows from tracelessness and orthonormality of matrices Tj. Comparing this 
relation with ([5]) and ([6]) we can write 

N 

Therefore, the length of the Bloch vector is not arbitrary. Bloch vectors V G M^^~^ lie on 
or within a sphere with radius Rl- Note that the space of Bloch vectors is different from 
the space of density matrices, although relation ([9]) establishes a one-to-one mapping 
between these spaces. 

3. Stratification of the Bloch space 

Boundary states are specified as such density matrices which possess at least one 
eigenvalue equal to zero. Let Pp denote a boundary state with p zero eigenvalues and 
with q = N — p nonzero positive ones: Ai,...,Ag (obviously Yl'j=i^j = 1) due to 
normalization). For such a boundary state relation ([5]) gives 

9 1 

i=i 

Due to lemma ffToj) we know that X]j=i -^j — ^/^- Therefore 
i2t 1 1 N — q p 

This inequality allows us to draw several conclusions. Some of them are obvious and 
well known, but we quote them just for completeness of the reasoning. Some other seem 
to be new. 

Boundary states with just one zero eigenvalue {p = 1) lie on or outside the small 
sphere (centered at pmax) of radius Rs = \Jl/N{N — 1). In other words, the small 
sphere contains density matrices with all N eigenvalues greater than zero. It is worth 
noting, that nonpositive but hermitian and normalized matrices can also be found just 
outside this sphere [8]. 

Pure states (with Ai = 1, all other eigenvalues are zeroes, so that p = N — 1) satisfy 
2 N-l 

d yPpurei Pmax) ^ "J^ • (-^^) 

Obviously, this result together with relation ^ imply that pure states must lie strictly 
on the surface of a large sphere [8j , so inequality becomes an equality. This explains 
why pure states are also extremal. Pure states must lie on the surface of the large sphere 
(which is a very well-known result) but the converse is not true. There are (on the large 
sphere) matrices which are not positive definite [8]. 

It seems to be a new conclusion that relation ( fT3l) implies that the space of density 
matrices is stratified by concentric spheres Sp (centered at Pmax) of radii 



p_ 

N{N-p)' 



rp=J ,rr.r p=l,2,...,iV-l, (15) 
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and boundary states Pp (with p zero eigenvalues) lie on or outside these spheres. We 
also note that for large and small spheres, we have Rs = ri and Rl = ttv-i- 

Lemma ( fT6j) allows us also to say that inequality (fT3|) is minimized for the boundary 
state of the form 

Pp = diagf-, -,0,...,o) = R{q), (16) 



q=N-p terms p terms 

SO only boundary states R{q = N — p) lie on the spheres Sp. Other boundary states 
with p zero eigenvalues lie outside these spheres. 

The discussed properties of the space of density matrices can easily be translated 
into the language of Bloch vectors. Firstly, density matrices with all nonzero eigenvalues 
correspond to Bloch vectors with the length 

IVP < ^ . (17) 

' ' N{N-iy ^ ^ 

Secondly, boundary states with p zero eigenvalues correspond to Bloch vectors such that 

|V|2 > 



,2 ^ P 



N{N -p)' ^ ' 

For boundary states R{q) defined in Eq. (fT6!) the above inequality becomes an equality 
and we have 

|Vh,„P = = ^ = -^jJ^. (19) 

All these lengths increase with growing number p. Finally, pure states correspond to 
|Vp = (A^ — 1)/A^ - on the large sphere. 

We have shown that the space of density matrices has stratified structure. With 
the growth of the number of zero eigenvalues the minimal distance of the given density 
matrix from the maximally mixed state also increases. 

4. Directional representation 

In any vector space, a vector can be represented by its direction and length. This is 
done writing 

V = |V|n, (20) 

with directional vector n being normalized: XljLi"^ — 1- Then, instead of represen- 
tation (j9]) we can take 

p = + |V|n ■ f = li^ + |V|Tfi, (21) 

where Tg = ^^^"^niTj. Representation (pTj) may be called a directional (or spherical) 
one. It was introduced and investigated by Kimura and Kossakowski [7]. Let us note 
that matrix T^i is hermitian and has the following properties 

Tr{Tfi} = 0, Tr{Tj} = 1, (22) 
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which follow from the tracelessness and orthonormality of Tj matrices. From the above 
relations it follows that matrix Tg possesses positive and negative eigenvalues which we 
will denote by {^k}k=i- They can be ordered in a natural way 

< yUl > /i2 > > /ijv < 0. (23) 

Subsequently, we can take the diagonal representation of the matrix in Eq. (!2T!) . The 
corresponding density matrix becomes also diagonal and its matrix elements become 

P., = ^ + |V|/i, J = 1,2,..., AT, (24) 

which must be nonnegative. The most restrictive relation is obtained for the smallest 
(negative) eigenvalue fiN = — |/^7v| (according to ( l23i) ). Inserting this eigenvalue into 
relation one finds that 

This estimate for the length of the generalized Bloch vector constitutes one of the main 
results of the work by Kimura and Kossakowski [7]. This is obtained here in a manner 
which seems to be simpler and more intuitive. Next, these authors investigated the 
minimal and maximal possible values of fit- They have obtained two cases 



'A^ - 1 1 
/ii = l^max = \/ jv ' ~~ ^{N 1) ' ^ = 2,3, . . . , A^. (26) 



/^fc = -y===, = 1,2, . . . , A^ - 1 and /^jv = l^min = ~y j^j- ^ ■ (^7) 



In the first case fl26|) . from inequality fl25|) we obtain 



N -1 

|V|<\/— (28) 
which is not unexpected. Next, from relations (l2l|) and (126|) we get 



'A^- 1 



1 






|V 


]V + 




1 






|V 


N ~ 





P.. = ^-|V|,/ ^^^^_^^ , ik = 2,...,N). (29) 

Allowing for maximum value of |V| as implied by (125]) we have 

Pu = l, Pkk = 0, {k = 2,...,N), (30) 

which corresponds to a pure state. 

The second case described by eigenvalues ( l27Il yields 
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so it corresponds to density matrices on or within a small sphere. Then, from relations 
(IM]) we find 



Pkk = ^r + \^\\ TTTTr TT' {k=l,...,N -1] 




P^^ = ^-|V|y^. (32) 

Taking maximum value of |V| from inequality (1311) we have now 

Pkk= (/c = 1, . . . ,iV - 1), p^^ = 0. (33) 

So we arrive at a boundary state with one zero eigenvalue lying on the surface of the 
small sphere. This summarizes the basic results of Kimura and Kossakowski [7J. On 
the other hand one may ask what interesting facts can be inferred. The answer is in 
the antipodean properties. This will be discussed in the next section. But before doing 
this, let us make two additional remarks. 

Firstly, we note that when inequality fl25|) becomes an equality then the correspon- 
ding density matrix ( pTI) is a boundary state. Indeed, in such a case, relation ( !24l) for 
j = N gives pnn = 0. Moreover, eigenvalue /ijv may be degenerate (see, for example 

m)- 

Secondly, let us construct a directional matrix Ta which describes the boundary 
state R{q) (Eq. fll6l) ). According to representation fl2T|) we now write 

Riq) = j^iN + \Vni,)\Tn. (34) 
Taking the length of Bloch vector as in f|T9|) and by definition f|T6l) we get 



V ^ 

q terms p=N-q terms 



as a directional matrix corresponding to the boundary state R{q). Note that the first q 
elements of this matrix are equal to |V/j(q)| as given in f|T9|) . 



5. Antipodean properties 

The antipodes are defined as follows. If a unit vector n defines certain direction, 
then vector (— n) specifies the opposite - an antipodean direction. Therefore, if a 
density matrix is given by its directional representation as in ( 1211) . then its antipodean 
counterpart is of the form 

Pa = ^In + |V|T_fi = - |V|Tfi, (36) 

because it seems obvious that T_fl = — T^. 
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Let us consider a state Pa{q) which is, by definition, antipodean to the boundary 
state R{q). Thus, we can write 



Pa{q) = ^i^v 



IV IT- 

I * a I -'■ n; 



(37) 



where is given in (135|) and |Va| is a suitably chosen length of the Bloch vector. Thus, 
matrix Pa{q) is of the form 



Pa{q) = diag 



N 




N-q 



qN 



N 



+ |V„ 



N{N-q)' 



{31 



q terms 



p=N-q terms 

The requirement of positivity entails a restriction on the allowed lengths of the Bloch 
vector. This implies 



IVJ < 



q 



N{N-qy 
Taking the maximal allowed value of 



(39) 



i?,(g) = diag(0,...,0,^ 



from (138!) we get 
1 



q 



N-q 



R{p). 



(40) 



g terms p=N-q terms 
which is antipodean with respect to the boundary state R{q). We note that any matrix 



satisfying inequality fl39|) is antipodean to R{q), while 
is the one most distant from the center, that is from the 



of the form of fl38l) . with V 
Ra{q) = R{p) given by Eq. go]) 
maximally mixed state pmax- 

This is a generalization of the conclusions derived by Kimura and Kossakowski [7] 
from the sets (!26l) and (!27|) of the eigenvalues of directional matrices. These two sets 
(up to numbering) clearly correspond to antipodal directions. In one direction there is 
a pure state (p = A^ — lorg = l) while in the opposite one we find at most the state 
(1331) with p=loTq = N-l. 



6. Final remarks 



We have shown that in the space of x A^-dimensional density matrices there are 
concentric spheres Sp (centered at the maximally mixed state) with increasing radii 
Tp = \Jp/N{N — p). These spheres characterize the stratification of boundary states 
Pp with p zero eigenvalues. States pp lie on the surface Sp or outside them. We have 
also shown, that only special boundary states R{q = N — p) defined in (fT6l) are on Sp. 
When a density matrix possesses q = N — p nonzero and unequal eigenvalues it must 
be outside the sphere Sp. The states situated inside spheres Sp have less than p zero 
eigenvalues, perhaps even none (all nonzero). In particular, p = 1 corresponds to a small 
sphere which contains only those density matrices which have all nonzero eigenvalues. 
Conversely, p = N — 1 correspond to pure states which lie only on the surface of the 
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large sphere. This shows that the space of density matrices has a complex stratified 
structure. 

The complexity of this structure is amplified by antipodean properties. If we assume 
that in a certain direction n we have a boundary state R{q), then in the antipode the 
corresponding density matrix Pa{q) is given by expression (|38l) . As long as the Bloch 
vector Va has the length satisfying the sharp inequality fl39|) we may expect that Pa{q) 
has all nonzero eigenvalues. When fl39l) becomes an equality, the antipodal matrix Pa{q) 
attains the form of a boundary state R{p = N — q) as given by Eq.( l40i) . So the point 
R{q = N — p) on the surface Sp corresponds to the antipodean point R{p = N — q) on 
the surface Sq. 

As yet, we are not fully aware what are potential consequences of the presented 
stratification of the space of N x N dimensional density matrices. Nevertheless, we feel 
that our findings are interesting and deserving further investigations. 
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